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Abstract 

We investigate the fundamental capacity limits of space-time journeys of information in mobile and Delay 
Q\ [ Tolerant Networks (DTNs), where information is either transmitted or carried by mobile nodes, using store-carry- 

forward routing. We define the capacity of a journey (i.e., a path in space and time, from a source to a destination) 
as the maximum amount of data that can be transferred from the source to the destination in the given journey. 
Combining a stochastic model (conveying all possible journeys) and an analysis of the durations of the nodes' 
O | encounters, we study the properties of journeys that maximize the space-time information propagation capacity, 

Qin bit-meters per second. More specifically, we provide theoretical lower and upper bounds on the information 
propagation speed, as a function of the journey capacity. In the particular case of random way-point-like models 
t — ■ (i.e., when nodes move for a distance of the order of the network domain size before changing direction), we show 

that, for relatively large journey capacities, the information propagation speed is of the same order as the mobile 
node speed. This implies that, surprisingly, in sparse but large-scale mobile DTNs, the space-time information 
propagation capacity in bit-meters per second remains proportional to the mobile node speed and to the size of the 
transported data bundles, when the bundles are relatively large. We also verify that all our analytical bounds are 
c/5 ' accurate in several simulation scenarios. 

O ' 

I. Introduction 

The problem of determining fundamental limits on the performance of mobile and ad hoc networks 
continues to attract the interest of researchers. Several important results have been achieved with the 
\ seminal papers by Gupta and Kumar [7] (which provided the first capacity bounds in static wireless 
networks) and by Grossglauser and Tse [6] (which showed that the mobility can increase the capacity of 
an ad hoc network). Various mobility models have been studied in the literature, and the delay-capacity 
relationships under those models have been characterized (e.g., [4], [13], [15]). However, the nature of 
these trade-offs is strongly influenced by the choice of the mobility model [14]. 

Moreover, there has been an increased interest in mobile ad hoc networks where end-to-end multi-hop 
paths may not exist and communication routes may only be available through time and mobility; depending 
on the context, these networks are now commonly referred as Intermittently Connected Networks (ICNs) 
or Delay Tolerant Networks (DTNs). Although limited, the understanding of the fundamental properties 
of such networks is steadily increasing. There is a significant number of results focusing on characterizing 
the packet propagation delay [3], [5], [17], assuming that packet transmissions are instantaneous, and more 
recently, the information propagation speed [8], [10], [11]. The authors of [3] took a graph-theoretical 
approach in order to upper bound the time it takes for disconnected mobile networks to become connected 

Part of this work will be presented in "On Space-Time Capacity Limits in Mobile and Delay Tolerant Networks", P. Jacquet, B. Mans 
and G. Rodolakis, IEEE Infocom, 2010. 

P. Jacquet is with INRIA, 78153 Le Chesnay, France. E-mail: philippe.jacquet@inria.fr 

B. Mans and G. Rodolakis are with Macquarie University, 2109 NSW, Australia. E-mails: bernard.mans@mq.edu.au, 
georgios.rodolakis@mq.edu.au 



X 



2 



through the mobility of the nodes. The papers [5], [17] analyze the delay of common routing schemes, 
such as epidemic routing, under the assumption that the inter-meeting time between pairs of nodes 
follows an exponential distribution. However, this assumption is not generally verified, depending on the 
relationship between the size of the network domain and the relevant time- scale of the network scenario 
under consideration [1], and this can result in either an over-estimation or an under-estimation of the 
actual system performance [2]. Departing from the exponential inter-meeting time hypothesis, in [10], 
[11], Kong and Yeh studied the information dissemination latency in large wireless and mobile networks, 
in constrained i.i.d. mobility and Brownian motion models. They showed that, when the network is not 
percolated, the latency scales linearly with the Euclidean distance between the sender and the receiver. The 
first analytical estimates of the constant upper bounds on the speed at which information can propagate 
in DTNs, again without considering the quantity of information that can be transmitted, were obtained 
in [8]. 

In contrast, in this paper, we investigate the space-time capacity of such networks, i.e., the maximum 
amount of information that can be transferred from a source to a destination over time. As the network is 
almost surely disconnected, we refer to journeys rather than paths, where a journey is an alternation of data 
transmissions and carriages using store-carry-forward routing. Informally, our objective is to determine 
how fast a given amount of data y can reach its destination. Formally, we use a probabilistic model of 
space-time journeys of packets of information in DTNs (in Section HI]), and define the journey capacity as 
well as the information propagation speed (in Section HIT]) , to provide the following main contributions: 

• we characterize the duration of node meetings, by bounding the probability function of the durations 
of the nodes' encounters, in Section HVt 

• we prove the first non trivial lower bounds on the information propagation speed (Theorem [[]), for a 
bounded journey capacity, in random waypoint-like mobility, in Section |Vl 

• we prove general upper bounds on the information propagation speed (Theorem [2] and Corollaries 1 
and 2), as a function of the journey capacity, and we investigate the properties of journeys that 
maximize the space-time network capacity in bit-meters per second, in Section |VIl 

• we compare and verify the analytical bounds with simulation measurements in Section IVIII 
We provide concluding remarks in Section IVIIIl 

II. Network and Mobility Model 

We consider a network of n nodes in a square area of size A = L x L and radio range R. As we want 
to focus on DTNs that are almost surely disconnected, we will analyze the case where R is fixed, while 
n, A — » oo, such that the node density v = ^ is bounded by some constant. 

Formally, we adopt the random geometric graph model [16]: two nodes at distance smaller than a 
maximum radio range R can exchange information. Moreover, we consider that the rate at which nodes 
can transmit data when they are within range is fixed, and equals G units of data per second. 
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Initially, the nodes are distributed uniformly at random. Every node follows an i.i.d. random trajectory, 
reflected on the borders of the square (like billiard balls). The nodes change direction at Poisson rate r and 
keep a constant speed v between direction changes. The motion direction angles are uniformly distributed 
in [0, 2ir) and are mutually independent among all nodes. When r = 0, we have a pure billiard model 
(nodes only change direction at the border). When r > 0, we have a random walk model; when r — > oo 
we are on the Brownian limit. When t — 0{j) — > we are on a random way-point-like model, since 
nodes travel a distance of order L before changing direction or hitting the border. 

III. Space-Time Journey Analysis 

We study journeys with a given capacity, i.e., journeys that guarantee that at least an amount of data can 
be transferred to the destination. Our aim is to find the shortest journey (in time) with journey capacity at 
least y, that connects any source to any destination in the network domain, in order to derive the overall 
information propagation speed. 

We base our analysis on a probabilistic model of journeys of packets of information that encapsulates 
all possible shortest journeys originating at the source, as used in [8]. Let C be a simple journey {i.e., a 
journey not returning to the same node twice). Let Z(C) be the terminal point. Let T(C) be the time at 
which the journey terminates. Let p(C) be the probability of the journey C. 

Let ( be an inverse space vector, i.e., with components expressed in inverse distance units. Let 9 be a 
scalar in inverse time units. We denote by w((,6) the journey Laplace transform, defined for a domain 
definition for ((,6): 

w(C,9) = E(eM-C-Z(C)-9T(C))) 

= £ c p(C)exp(-C-Z(C)-0T(C)). 

We call p(z ,z 1 ,t) the normalized density of journeys starting from z at time 0, and arriving at z 1 
before time t: 

p(z ,z u t) = — P( C ) ■ 

\\z!-Z(C)\\<R,T(C)<t 

Let us consider that a bundle of information of y bits is generated at t = on a node at coordinate z = 
(x ,y ). Let us initially consider a destination node which stays motionless at coordinate z 1 = (xi,yi); 
in this case, p(z , z 1 , t) denotes the probability that the destination receives one bit of information before 
time t. Now, let us consider a moving destination node, that is located at coordinate z 1 = (x 1: yi), at time 
t. We denote z = z\ — z . Let q(z,t,y) denote the probability that there exists a journey of capacity at 
least y reaching the destination before time t. 

The information propagation speed s(y), considering a journey capacity y, is defined as the minimum 
ratio of distance over time above which the journey probability tends to 0, i.e., 

. if UfU > s(y), then \im M ^ 00 q(z,t,y) = 0; 
# if < s ^ then li m || z || t ^ oo g( Zj t ? y) > o. 

We also define the space-time information propagation capacity c(y) (from now on simply referred 
to as the space-time capacity), as the maximal transport capacity in bit-meters per second, that can be 
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achieved by any journey of capacity y. Thus, in this model, the space-time capacity corresponds to the 
product c(y) = s(y)y. 

Therefore, in order to determine the space-time capacity limits of mobile and delay tolerant networks, 
we will analyze the information propagation speed, as a function of the journey capacity; in the following 
sections, we will compute lower and upper bounds. In order to derive the bounds, we first study the 
characteristics of node meetings. 

IV. Node Meetings 

A meeting (or encounter) between two nodes occurs when their distance becomes smaller than or equal 
to R, i.e., when the nodes come into communication range. 

Lemma 1: A node A, moving in direction ip , meets new nodes moving in direction between ipi and 
ipi + dip at rate: \ ^ = sin ^ 1 "^ )dijj, for ip , ipi G (— 7r, 7r], where R is the radio range. 

Proof: See appendix. ■ 

We denote the meeting duration by the random variable T. 

Lemma 2: The probability P(T > t) that a meeting has duration at least t satisfies: 

P(T >t)<min(l,^). 

8vt 



Proof: The average number of neighbors of any node is ixvR 2 . From Lemma [Q the rate at which a 
node meets new neighbors is / = Therefore, from the Little formula, the average meeting time {i.e., 
the time that a node remains a neighbor) equals ZL y- = The proof follows by applying Markov's 
inequality. ■ 
In the pure billiard model (i.e., when r = 0), we can give the exact formulas on the meeting time 
distribution. We note that our model where nodes bounce on the borders like billiard balls is equivalent 
to considering an infinite area made of mirror images of the original network domain square: a mobile 
node moves in the original square while its mirror images move in the mirror squares [8]. 

Lemma 3: We denote the meeting duration by the random variable T. The probability density function 
pr{t) of T is: 



Pr(t) = | log 



¥ + 1 



R 2 \ R 



for t > 0, where v is the node speed, R is the radio range. 
When t — > oo, the cumulative probability P(T > t) is: 

P(T>t) = -P-, + 0' '' 



3(vt) 2 \(vt) 



Proof: See appendix. 
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V. Lower Bound 



We prove a lower bound s^(j/) on the information propagation speed, for journey capacity y, in the 
random way-point-like mobility model, i.e., when nodes travel a distance of the order of the network 
domain length before changing direction. Initially, we focus on the pure billiard mobility model, i.e., we 
assume that nodes do not change direction unless they hit the border. Finally, we remark that the result 
can be generalized to node mobility with a small change of direction rate. 

We will show that, for all destination nodes which, at time t, are at distance r ~ s L (y)t of the initial 
source location, there is a journey of duration t and of capacity y from the source to the destination, 
with probability strictly larger than 0. We consider large distances r = <d(y/n), where n is the number 
of nodes in the network; in this case, the square network domain has a side length r = 0(-y/n), as we 
are interested in the case where the node density is constant (but strictly larger than 0), as discussed in 
Section lU We show that, when the journey capacity is y < — , for a constant K, the lower bound is 
s l(u) = v, where v is the mobile node speed. 



Fig. 1. Definitions of rendez-vous point A of the information generated at location S with the destination D (left), and of angle (j>c with 
respect to the speed of node C and location B (right). 

We consider a source node S and a destination node D. We denote by v s and v D the respective vector 
speeds of the source and the destination. We assume that the source starts sending the information at time 
0. We define the point A as the third vertex of the isosceles triangle, formed with the two other vertices 
located at S and D (at time 0) and with sides SA and DA of equal length r, while DA is parallel to 
the destination speed v^, as illustrated in Figure [Q Point A is therefore the rendez-vous point of a node 
moving at constant speed v, in the direction of SA, and the destination node, while the nodes contact (at 
the same location) occurs at time t A — -. Similarly, if the (asymptotic) information propagation speed 
is equal to the node speed v, the information will reach the destination at location A' = A ± AZ, with 
\AZ\ = o(r), at time t A > =t A + o(~). 

We will describe a routing scheme that constructs a journey of duration t A = r + o(-), which originates 
at S and ends at any given point A, and guarantees that for any direction of the destination node speed, 
the journey capacity is at least y. We assume w.l.o.g. that the radio range is R = 1 and the communication 
rate is also G = 1, to simplify the expressions (to generalize, it is sufficient to perform a simple scaling). 
We note that, in this case, ensuring a journey capacity at least equal to y is equivalent to ensuring a 
minimum meeting duration y for all transmissions in the journey. 

The routing scheme proceeds in three stages, illustrated in Figure [2l In all stages, the information is 
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Stage 1 




Fig. 2. Overview of the routing scheme achieving the lower bound of information propagation towards the rendez-vous point A, in three 
stages. 



passed among nodes moving at relative direction of angle between | and a, with a value of a that we 
will precise in the following. 

Initially, we consider a point B located on the destination's trajectory (before the rendez-vous point A). 
We also take B such that the distance from the rendez-vous point A is r B = 0(- v /r). In the first stage, the 
information is transmitted to new nodes (according to the above angle restriction and ensuring a journey 
capacity at least y) until reaching a node, whose trajectory's distance from B is at most y/r. 

In the second stage the node with the information simply travels a straight line (of length r + 0(\/r)) 
until approaching the point B within distance y/r. 

In the third stage, the information is transmitted to new nodes (again, with a relative direction angle 
in [| ,a), and ensuring a journey capacity at least y) until the information is transmitted to a node that 
passes within distance 1 of the rendez-vous point A, while the contact duration with the destination is 
sufficient to transfer all the information. 

We will show that this routing scheme guarantees that the information will reach the destination with a 
journey of capacity at least y, with a total journey duration of - + 0(— ). More precisely, we show that 
the duration of the first and third stages is 0(— ). Since the duration of the second stage is L + 0(— ), 
a lower bound on the information propagation speed is v. 

We now analyze the duration of the three routing stages. 

1) Stage 1: We introduce the following notations. Let C be the node that most recently received all 
the information, moving at speed v c . We define <pc as the angle formed between the vector CB (defined 
by the locations of the node C and the point B) and the speed v c , as depicted in Figure [Q 

Lemma 4: The duration t\ of stage 1 of the routing scheme is 6(— ), almost surely. The distance 
traveled is 0(y/r). 

Proof: See appendix. ■ 

2) Stage 2: 

Lemma 5: The duration t 2 of stage 2 of the routing scheme is - + 0(— ), almost surely. 

Proof: The initial distance SB is at most r + r B = r + 0(\/r). From LemmalU the distance T\ = CA 
at the end of stage 1 is r + 0(y/r). The minimum distance of node C trajectory to B, and is at most 
r 2 = ri sin(4^) = y/r + 0(r~^), as depicted in Figure [2l Therefore, there is a point in the trajectory such 
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that the final distance of node C from the point B is exactly y/r. Therefore, the total distance traveled in 
stage 2 is at most n(l + (^)) = r + 0(y/f). ■ 

3) Stage 3: At the beginning of stage 3, there is a node carrying the information, located within distance 
v b + \pf' from the rendez-vous point, and within distance yfr from the destination's trajectory. In this 
stage, the information is transmitted to new nodes (again, according to the above angle restriction and 
ensuring a capacity at least y) until reaching a node that passes within distance 1 of the rendez-vous point 
A, while the contact duration with the destination is at least y. 

Equivalently to stage 1, let C be the node that most recently received all the information, moving at 
speed vc- We introduce again the angle 4>c, this time defined with respect to the rendez-vous point A; 
namely, <pc is the angle formed between the vector CA (defined by the locations of the node C and the 
rendez-vous point A) and the speed v c . 

Lemma 6: We consider a node C, at distance rc from the rendez-vous point, moving with speed 
at a direction such that the relative angle with the destination's direction is at most a = If the angle 
(j) C is at most then the trajectory of C passes within range of the destination and guarantees that the 
meeting duration with a destination located at A, moving at constant speed, will be at least equal to y. 

Proof: The relative speed of the node C, with respect to the destination's speed, is at most 2v sin(|) < 
va. If the node C passes within distance m from the rendez-vous point, the meeting duration is at least 
(since the distance traveled within range, in the frame of reference of the destination, is at least 
1 — m). Therefore, in order for the meeting duration T to be at least equal to y, it is sufficient that: 
m < 1 — yva — \. In this case, we guarantee a meeting duration at least equal to y. Moreover, if we have 
4>c < me n °d e W ^ P ass w i mm distance | from the rendez-vous point. ■ 

Lemma 7: The duration t 3 of stage 3 is O(^), almost surely. At the end of stage 3, the destination is 
reached at the rendez-vous point with probability strictly larger than 0. 

Proof: See appendix. ■ 

Theorem 1: Consider a network with constant node density v, radio range R and communication rate 
G, where nodes move at speed v > and change direction at rate r = 0. When the journey capacity is at 
most y — — , where K is a constant, a lower bound on the information propagation speed is s L (y) = v. 

Proof: Considering the final position of any destination, we can define a rendez-vous point A. If 
the distance of the rendez-vous point from the source location at time is r — > oo, based on the previous 
lemmas, there exists with strictly positive probability a journey of capacity at least y that reaches any 
rendez-vous point A within time ~ L . Therefore, the asymptotic information speed is at least v. ■ 

We note that, in case the network domain A = L x L is sufficiently large, for all destination nodes 
which, at time t = Q(L), are at distance r = o(vt) of the initial source location, there is almost surely a 
journey of duration t and of capacity y from the source to the destination. 

Remark 1: Although, we derived the lower bound in a pure billiard mobility model, the proof can be 
easily generalized to a random walk model, where the change of direction rate is O(^), by restarting from 
the first stage at any change of direction (an event which occurs a finite number of times). 
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VI. Upper Bound and Space-Time Capacity 



In this section, our aim is to find the shortest journey of capacity at least y that connects any source to 

any destination in the network domain. We prove an upper bound su(y) on the information propagation 

speed, for journeys of capacity y. 

Theorem 2: Consider a network with n mobile nodes with radio range R, communication rate G, in 
a square area of size A = L x L, where nodes move at speed v, and change direction at rate r. When 
n —y oo, such that the node density becomes v — -p, an upper bound on the information propagation 

speed, for journeys of capacity y, is the smallest ratio of - with: 



Remark 2: The expression of 9 has meaning when ixvR 2 ^(y) < 1. Above this threshold, the upper 
bound for the information propagation speed is infinite. Such a behavior is expected, since there exists 
a critical node density above which the graph is fully connected or at least percolates [12]. In addition, 
according to Theorem [2l in percolated networks, there is a critical journey capacity y c , such that, when 
y > y c , the propagation speed is bounded by a constant. 

Proof: We assume that a source starts emitting information at position z = and time t — 0. We 
consider the probabilistic space-time journey model presented in Section (Till which includes all shortest 
journeys originating at the source. Equivalently, we model journeys of very small beacons of information, 
such that beacon transmissions are instantaneous. 

We initially consider an infinite network with a Poisson density of nodes A. We will upper bound 
the probability density of journeys in the infinite network model. However, by applying an analytical 
depoissonization technique [9], we obtain an equivalent asymptotic estimate of the journey density when 
the number of nodes n is large but not random. 

We decompose the journeys into two types of segments, modeling node movements and beacon trans- 
missions: 

• emission segments S e (u, v): the node transmits immediately after receiving the beacon; v is the speed 
of the node that just received the beacon, and u is the emission space vector and is such that |u| < R; 

• move-and-emit segments S m (u, v, w) = M(v, w) + u: M(v, w) is the space-time vector correspond- 
ing to the motion of the node carrying the beacon, where v is the initial vector speed of the node 
when it receives the beacon and w is the final speed of the node just before transmitting the beacon; 
the vector u is the emission space vector which ends the segment. 

Considering any sequence of segments, we can always upper bound the segment probabilities (see [8], 
Section III-B). In fact, the conditional probabilities, given the node direction and speed, are upper bounded 
by unconditional probabilities: 
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• P(S e (u)) = P(u)A, where P(u) is the probability density of u inside the disk of radius R, and A 
is the node density (to make the emission possible); 

• P(S m (u, v, w)) = P(u)P(M(v, w))2v\, where P(u) is the probability density of u on the circle 
of radius R (we only need to consider the earliest transmissions, which occur at the maximum radio 
range), P(M(v, w)) is the probability that the node movement equals the space vector M(v, w), 
and v is the node speed. 

This upper bound journey model results in a higher density of journeys than in the actual network. But, 
in this model, any journey can be decomposed into a sequence of independent segments. Consequently, 
we can express a journey C as an arbitrary sequence of emission or move-and-emit segments, i.e., using 
regular expression notation, C = (S e + S m )*. Moreover, we can calculate the Laplace transform of the 
journey probability density, based on the Laplace transforms of the segments. We denote the segment 
Laplace transforms by l e ((,9) = E(e~^ ,9 ^' Se ) and l m ((,9) = E{e"^ ,e " ) ' Sm ), for emission and move-and- 
emit segments, respectively. Equivalently to the formal identity = l+x+x 2 +x 3 + ..., which represents 
the Laplace transform of an arbitrary sequence of random variables with Laplace transform x, the journey 
Laplace transform has a denominator k((,9), equal to: 

fc(C,0) = i-(k(C,0)+UC,0))- (2) 

We have the following Laplace transform expressions: 

• l e {Qi 9) = E(e~^' u ), where u is uniform in the disk of radius R, with density A, i.e., l e (£,9) = 
AvrffJi(KlP). 

. l e ((,9) = P(e-^ u )P(e~ (<: ' e) ' M(v ' w) ), where u is uniform on the circle of radius R, with density A, 
i.e., E(e-^) = 2vrAPJ (|C|P), and £( e —^(v,w)) = 1 (see [8]). 

We derive an upper bound on the information propagation speed, in the special case where the journey 
capacity is y = 0, from the analysis of the singularities of the journey Laplace transform, for A equal to 
the node density in the network (cf. Theorem 1 in [8]). The upper bound is the smallest ratio - of the 
non-negative pair (p, 9) which is a root of the denominator k(p, 9) (with p = \(\), obtained by substituting 
the segment Laplace transforms expressions in ([2]). 

In order to generalize to journeys of a given capacity y > 0, we will restrict the set of possible 
journeys, to those satisfying the desired capacity constraint, and calculate the Laplace transform of the 
journey density in this restricted set. 

First, we remark that, a journey has a capacity at least y, if and only if the journey thickness (i.e., the 
minimum duration of all data transmissions in the journey) is at least equal to ^, with G the communication 
rate. Therefore, when considering journeys of a given capacity, we can equivalently focus on the possible 
journeys with minimum node meeting duration ^. 

Therefore, in the upper-bound journey model, we can substitute the probability of any emission segment 
with the probability of the same emission segment, while additionally ensuring that the emission duration 
is at least ^. Thus, for the singularity analysis, we substitute in (O the Poisson density A with a node 
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density v^{y), where j(y) is an upper bound on the probability that the meeting duration is at least ^. 
This direct substitution is feasible because we work with the upper bound journey model, where successive 
segments (including all transmissions) are independent of the previous network state. Again, this results in 
considering a higher density of journeys than in the actual network (including some transmissions which 
are not actually possible, due to the node directions); however, there is no impact on the validity of our 
analysis, since we are interested in upper bounds. 

To conclude the proof, it suffices to substitute quantity 7(1/) using Lemma |3] when r = 0, and Lemma |2] 
when r > 0. ■ 

We derive the following corollaries expressing the behavior of the upper bound when the journey 
capacity is large, in random waypoint-like (r — > 0) and random walk/Brownian motion mobility (r > 0), 
respectively. 

Corollary 1: When nodes move at speed v > 0, and ^ — > (i.e., the journey capacity y is large) with 
r > 0, the propagation speed upper bound is O(^J^-Rv). 

Proof: See appendix. ■ 
Corollary 2: When the node speed is v > 0, and ^ — ► (i.e., the journey capacity y is large) with 
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r = O(j) — > 0, the propagation speed upper bound is v + 0(r + 

Proof: See appendix. ■ 
We observe that, for large journey capacities y, the upper bound on the information propagation speed 
su(y) tends to the actual mobile node speed v in random way-point-like mobility, while it decreases with 
the inverse square root of the journey capacity y in random walk or Brownian motion mobility. In both 
cases, the resulting upper bound on the space-time capacity c(y) = s v (y)y is a function which increases 
with y. 

Remark 3: When nodes move at speed v > in random way-point-like mobility: 

• from Theorem [Q a lower bound on the propagation speed is v, for any bounded y, and when the 
node density is v = 0(1); 

• from Corollary [2l an upper bound on the propagation speed is v, for journey capacities y such that 

v = o(y). 

Therefore, we notice that there is a range of values of y, for which our bounds are almost tight. More 
generally, we deduce that the information propagation speed in random way-point-like mobility models 
is of the same order as the mobile node speed, for (bounded) journey capacities that are relatively large 
with respect to the node density. 

This implies that, in sparse but large-scale mobile DTNs, the space-time information propagation 
capacity in bit-meters per second remains proportional to the mobile node speed and to the size of the 
transported data bundles, when the bundles are relatively large. It is rather surprising that the propagation 
speed does not tend to when the size of the bundles increases, which would result in a sub-linear 
increase of the space-time capacity. 



11 




Fig. 3. Snapshots of simulated information propagation at three different times (t — 100, 170, 240), for a small journey capacity y = 0.5 
(top) and a larger journey capacity y — 2.5 (bottom). Larger black squares represent nodes that have received all the information at the time 
of the snapshot. 

VII. Numerical Results 

In this section, we perform simulation measurements to compare to the analytical bounds on the 
information propagation, derived in the previous sections. We developed a simulator that follows the 
network and mobility model described in Section HEl We simulate the epidemic broadcast of information, 
and we consider journeys with a given lower bound on the capacity y, as described in Section [TTTJ We note 
that the simulation is more general than the simple broadcast of a packet of size y, since the information 
can also be transferred on a given journey using smaller packets. In fact, we precisely ensure that the 
journeys of the simulated broadcast have a capacity at least y, without imposing further restrictions. For 
all the following simulations, we consider a communication rate G — 1 units of data per second (e.g., 
if one unit of data corresponds to x Mbits, the journey capacity in the following examples should be 
multiplied by x Mbits). 

We first show how information propagates in a full epidemic broadcast, by illustrating two typical 
and distinct situations, depending on the journey capacity y. In the simulated scenario, a source starts 
broadcasting information at time t = 50, in a network of 5000 nodes, in a 2000m x 2000m square, with 
radio range R = 10m, and mobile node speed v = 5m/ s, with pure billiard mobility (r = 0). In Figure [3] 
we consider two cases: a smaller journey capacity y = 0.5 (top) and a larger journey capacity y = 2.5 
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(bottom). For each case, we depict three snapshots of the simulated information propagation at three 
different times, t = 100, 170, 240, from left to right. The small black dots represent the mobile nodes; 
when two dots are in contact, the corresponding nodes are within communication range. The larger black 
squares represent nodes that have received all the information at the time of the snapshot, i.e., those that 
can be reached by a journey of capacity y. The simulation scenario is exactly the same in both the top 
and bottom figures, with the only change concerning the journey capacities. In both cases, the location 
of the source is approximately located at the center of the disk containing the black squares, at the top 
left figure. We observe that, at the top row of Figure [3] corresponding to a small journey capacity, the 
information propagates as a full disk that grows at a constant rate, which coincides with the information 
propagation speed; all nodes inside the disk can be reached by a journey of capacity y, almost surely. 
Equivalently, this means that the average information propagation delay scales linearly with the distance 
from the source, and the ratio of the propagation delay over the distance is equal to the inverse of the 
information propagation speed. On the other hand, at the bottom row, corresponding to a larger journey 
capacity, only some of the nodes inside the disk have been reached by a journey of capacity y. In this 
case, the average information propagation delay does not necessarily scale linearly with the distance from 
the source. However, the information still propagates at a (smaller than before) maximum speed, equal to 
the rate at which the disk radius grows. 

Next, we simulate a network of 500 nodes, moving in an area 600m x 600m, with a radio range of 10m, 
a mobile node speed of 5m/s and a communication rate G = 1 units of data per second. We simulate 
two different mobility parameters (rates of direction change): r = for the pure billiard mobility model, 
where nodes change direction only when they bounce on the border, and r = 0.05 for a random walk 
model. 

In Figure |4] we plot the ratio of the propagation delay over the distance from the source, versus the 
distance, for journey capacities y = {1; 2; 3}. Each sample point in the plots corresponds to a simulation 
measurement. The distance is measured from the location of the source when the information was emitted 
to the location of the destination when the information was received. We notice that, for all journey 
capacities, the ratio of the propagation delay over the distance is larger than a non-zero constant. The 
constant lower bound on the ratio, in this simulation scenario, is close to the inverse of the mobile 
node speed (which is plotted in the figures as a straight line, for comparison). Furthermore, this constant 
corresponds to the upper bound on the information propagation speed, which was calculated in Theorem [2] 
In fact, for small journey capacities (e.g., y = 0.5), we notice that the upper bound on the information 
propagation speed is larger than (but close to) the mobile node speed. For larger journey capacities and 
r = 0, the upper bound can be obtained from Corollary [2l and indeed corresponds to the mobile node 
speed. We also notice that, for r = 0.05, the average distance that each node travels before changing 
direction is 100m, which is of the order of the square network domain length. Therefore, in this case, 
the upper bound on the propagation speed also remains close to the estimate for random waypoint-like 
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Fig. 4. Ratio of information propagation delay over distance versus distance from the source, for different journey capacities (y — {1; 2; 3}, 
respectively), compared to the inverse of the mobile node speed, with pure billiard mobility (r = — left), and random walk mobility 
(r = 0.05 - right). 



mobility in Corollary [2l i.e., the mobile node speed. 

In Figure [51 we depict the simulated average propagation time versus the distance, for several different 
journey capacity values y = {0.5; 1; 1.5; 2; 2.5; 3}. Time is measured in seconds, and distance in meters, 
therefore, the inverse slope of the plots provides us with the information propagation speed in ms^ 1 . We 
compare it to a line of fixed slope corresponding to the mobile node speed. For comparison, we plot the 
theoretical upper bounds on the information propagation speed (derived from Theorem [2]) in Figure [6] 
Simulations show that the theoretical speed is clearly an upper bound. Moreover, we notice that the upper 
bound in the case corresponding to random waypoint-like mobility is tighter, due to the fact that our 
analysis of the node encounter duration analysis (see Lemma |2]> is exact in this case. 

In Figures [51 we also notice that, for journey capacities up to 2 units of data per second, the measure- 
ments rapidly converge to a straight line of fixed slope, which implies a fixed information propagation 
speed, as illustrated by the top row of Figure [H However, for larger journey capacities, border effects 
become significant and the slope of the measurements tends to 0; this means that, although the maximum 
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Fig. 5. Average propagation delay versus distance for different journey capacities (y — {0.5; 1; 1.5; 2; 2.5; 3}), with pure billiard mobility 
(r = — top), and random walk mobility (r = 0.05 — bottom). 
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Fig. 6. Upper bound for the information propagation speed as a function of the journey capacity (n = 500, A = 600m x 600m, R = 10m, 
v — 5m/ s, G = 1 units of data per second), with pure billiard mobility (r = — left), and random walk mobility (r = 0.05 — right). 



information propagation speed is still a non-zero constant, the information does not propagate uniformly as 
a disk growing at constant speed. In this case, information propagation occurs similarly to the expectation 
illustrated in the bottom row of Figure [3j 

Finally, in Figure |7J we plot the space-time capacity in bit-meters per second, versus the distance from 
the source, achieved by journeys of different capacities y = {0.5; 1; 1.5; 2; 2.5; 3}, in the same simulation 
scenario. The space-time capacity is obtained by multiplying the average propagation speed s(y) with 
the journey capacity y. We observe indeed that, for journey capacities up to 2 units of data, the plots of 
the space-time capacity in Figure [71 converge to c(y) = s(y)y ~ vy; this is consistent with Remark [3] 
For larger capacities, the space-time capacity has not converged to a constant value, due to the fact that 
the network domain is finite. However, we note that, in a larger network, the space-time capacity would 
be larger for journeys of larger capacities. In fact, in an infinite network, the space-time capacity would 
converge to a constant value for any finite journey capacity. 
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Fig. 7. Space-time capacity in bit-meters per second, versus distance from the source, for journey capacities y = {0.5; 1; 1.5; 2; 2.5; 3}, 
with pure billiard mobility (r = — top), and random walk mobility (r = 0.05 — bottom). 



VIII. Concluding Remarks 

We characterized the space-time capacity limits of mobile DTNs, by providing lower (Theorem [D and 
upper bounds (Theorem [2]) on the information propagation speed, with a given journey capacity. Moreover, 
we verified the accuracy of our bounds with extensive simulations in several scenarios. 

Such theoretical bounds are paramount in order to increase our understanding of the fundamental 
properties and performance limits of DTNs, as well as to design or optimize the performance of specific 
routing protocols. In fact, our results provide lower and upper bounds on the best achievable propagation 
delay of bundles of data, over large distances. 

It is also worth noting that our analysis provides the first known lower bounds on the information 
propagation speed in mobile DTNs (for random waypoint-like mobility models), and generalize previously 
known upper bounds. 

More specifically, in the case of random waypoint-like mobility models, we showed that for relatively 
large journey capacities, the information propagation speed is of the same order as the mobile node speed. 
This implies that, in sparse but large-scale mobile DTNs, the space-time information propagation capacity 
in bit-meters per second remains proportional to the mobile node speed and to the size of the transported 
data bundles, when the bundles are relatively large. 
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Appendix 

A. Proof of Lemma [7] ( Meeting Rate ) 

When n, A — > oo, we can consider an infinite network with a Poisson density of nodes v — to 
simplify the proof. In fact, if we consider an area A of the infinite network, the number of other nodes 
is given by a Poisson process of rate n, and we can depoissonize it [9], to obtain the equivalent result 
when the number of nodes n is large but not random. 

Let u be a unit vector always centered at the position of node A. We denote by / the rate at which 
mobile nodes enter the neighborhood range of node A at position Ru with respect to the node location 
ZA(t), where R is the radio range. 

Let us denote by B a second network node, with a constant vector speed \ B . The Poisson density of 
presence of B at any location on the plane is v. The relative speed of the nodes is \ B — V4. The projection 
of the relative speed on the vector _Ru equals (Ru ■ (vg — v^)) u. The rate at which any node B enters 
the neighborhood range of the node A at u, is /(v^, vb, u) = max{0, u ■ (vb — va)vR}- 

By averaging on u, we have the total meeting rate: 

\vb - va| cos^uRdip = 2v\v B - v A \R- 

■f 

Therefore, the rate at which a node meets new neighbors is proportional to their relative speed. From the 
law of sines, the relative speed is proportional to sin(4^), where Atp = ipi — tp Q is the angle formed 
between the speed vectors. By normalizing, we obtain the meeting rate. 

B. Proof of Lemma \3\ ( Distribution of Encounter Duration ) 

We consider the encounter of two nodes A and B, moving at speeds v^ and v B respectively. We define 
Av = vb — \a as the relative speed of the nodes. Therefore, taking as a frame of reference the position 
of node A, node B is moving at constant speed Av, as illustrated in Figure [T0l We denote by Av the 
Euclidean norm of the relative speed (i.e., the relative velocity). From the law of cosines, it holds: 

Av = \\v B - v A \\ = 2vsm(—), (3) 

where ip G [0, 2tt) is the angle between the node speed vectors. 
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From Lemma [TJ the rate at which nodes meet is proportional to their relative speed. Therefore, 
normalizing ©, the angle tp is distributed according to the probability density function: 



p^(x)dx = -sin(—)dx, x G [0, 27r), 



(4) 



and, substituting V = 2t>sin(|) according to ©, the density function p& v (V) of the relative velocity is: 



PAv(V)dV 



V 



V 2 



dV, V e [0, 2v\ 



(5) 




Fig. 8. Encounter of nodes A and B in the frame of reference centered at A: Av is the relative speed of B, d is the length of the chord 
traveled by B within range, £ is the distance of the chord d from A. 



Fig. 9. Probability density function pr(t) = | log ffrr (1 



(vt) 



■) — i of the node encounter duration T, for v = 1. 



Always in the frame of reference of node A, we denote by d the distance traveled by node B within 
range of node A. In other words, d is the length of a chord of the circle of radius R (the radio range), 
centered at node A. We define t as the distance of the chord from A, as depicted in Figure [TOl We remark 
that, as a node moves and meets new neighbors, quantity i is distributed uniformly at random between 
and R, since meetings occur equiprobably at any point of the diameter perpendicular to the node relative 
speed. Therefore, since d = 2V R? — £ 2 , the distribution of the length d is: 



P(d > x) 



1 - 



x z 

4R2' 



Differentiating, we obtain the probability density function: 



X 



2R^AR 2 - x 2 



x e [0,2R]. 



(6) 



(7) 
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If T is the duration of the encounter, we have: 

d = Av x T, 



(8) 



where all quantities are random variables. 

Let us consider a given relative velocity Av = V. In this case, we can define the conditional probability 
density px(t \ Av = V) of the encounter duration, with t E [0, y}: 

pr(t | Av = V)dt = p d (x | Av = V)dx = p d (Vt)Vdt, 



where x = Vt, according to ®. 
Combining with ©, 



pr(t | Av = V) 



VH 



2y/A-(Vty' 

Considering the probability density function pr(t), and using © and ©, we have for t > 0: 



(9) 



p T {t) 



n2v 

/ pr(t | Av = V) x PAv (V)dV 
Jo 



I lo s 



i + 



R 2 



vty 



R_ 

2vt' 



We note that the fact that nodes bounce on the borders does not impact on this result. We plot the 
probability density function pr(t) (for R = 1, v = 1) in Figure [9j 
By simple integration, we obtain the probability P(T > t): 



P(T>t) = - log 



¥ + 1 



R v 
vi~ R 



1 

2' 



(10) 



For large t, we have 4^ > 0. Therefore, using the identity log 2 = 2J2^ =0 -^zi (f^) Z ™ +1 ' we have: 



x -l\2n+l 



P{T > t) 



R2 +0( R ' 



3(vt) 



vty 



C. Proof of Lemma |?] (Duration of Routing Stage 1 ) 

Since we consider meetings of relative angle at most a, the relative speed of two meeting nodes is 
maximized when the angle between them is a (and equals 2d sin (|)). Therefore, in order for the meeting 
duration T to be at least equal to y, it is sufficient that the distance d traveled within range, in the frame 
of reference of one of the nodes (see Figure [TO]), satisfies: 



d > vay > 2vsm{—)y. 



According to ©, P(d > x) 



1 - f , and P(T > y) > Jl - Assuming that y > i we 



take a = 2^' 



P(T>V)>^>1 



For smaller y, the same bound clearly still holds. 
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Fig. 10. Encounter of nodes A and B in the frame of reference centered at A: Av is the relative speed of B, d is the length of the chord 
traveled by B within range, t is the distance of the chord d from A. 



From Lemma [H the probability to meet a node at an angle in [|,a] is P a = cos(|) — cos(|) > j^a 2 , 
since a < |. 

The rate at which a node meets new nodes at such an angle, ensuring that the meeting duration is at 
least y, is: 

/. > —KP(T > V) > 

7T 16 

We note that the angle cpc determines the distance 4 of node C trajectory from the point B (see 
FigureQ]). In fact, it holds: ds = \CB\ sin <p C - When a node moves, <p c varies, while ds remains unchanged. 
In fact <pc always increases when a node moves towards the destination. However, after a node movement 
of distance 5, we have A<fi c = O(j^), and if 5 = o(\CB\), 4> c is not modified asymptotically. 

Thus, if the initial angle between the source and the destination is b, the expected time E(t\) until 
f < <\>c < « is: 

„, , N 2b 32tt n/ 1 . 

E(t[) < — < — = e — . 

aj a 6 vv vv 

From Lemma [U the rate at which a node meets nodes at relative angle [ip,ip + dip] is — sin(|)c?/3. 
Therefore, the node C that last received the information meets new nodes C with angle (j)' c < and 
with meeting duration at least y, with rate (assuming that <\>c remains between | and a): 



/a > ^(T > y / sin - + x)dx > 1. 4> + O r"' 

V r 

and the expected time E(t") until meeting such a node is 0(— ) (we note that - < 2K). We notice that 
the t'l = o(r) almost surely, and we can indeed assume that 4> c remains constant until meeting C . 

Therefore, it holds that the duration t\ of stage 1 is t\ — t[ +t" = 0{—) almost surely. The distance 
traveled is vti+0(-), where the second term corresponds to the further distance moved by the information 
in O(-) transmissions. Since - = 0(1), the total distance traveled is 0(— ). 

D. Proof of Lemma ( Duration of Routing Stage 3 ) 

We proceed equivalently to stage 1. Stage 3 ends when a node with angle <\>c < ^ receives the 
information. Equivalently to the proof of Lemma HI the expected time t' 3 until the relative speed of the 
node to the rendez-vous point A is between | and a is E(t' 3 ) = 0(— ). 
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We consider meetings with nodes C, such that 2r c > < t^, where k\ > is a constant. The node C 
that last received the information meets new nodes C with angle <pc> < % (< 55^7)' an( ^ w ^ m meetm g 
duration at least y, with rate (assuming that 4>c is between | and a): 

/ 2 > P T >!/)/, sin - + arjdx > - u ' + O r"' 

V r 

and the expected time E(t'l) until meeting such a node is 0(^7). Since C varies, if it becomes larger 
than a (or smaller than |), the information is forwarded to a new node such that 4> c is between | and a 
again (in constant time). 

Moreover, we have indeed that r c = 0{^fr) < ^(1 + 0(1)) for some positive constant k\, since the 
distance traveled at this stage is at most vt' 3 + vt 3 + 0(-) = 0(y/r). We assume that r c < r B , which 
we can ensure by choosing point B sufficiently far from the rendez-vous point A. In this case, when 
the information is transmitted to node C, the node's direction, with respect to the destination's speed, is 
of angle at most 0c* < a. Therefore, after time t 3 = 6 (-^7), the destination is reached with probability 
strictly larger than 0. 

E. Proof of Corollary [7] 

W. 1. o. g., we take R = 1 and G = 1. Let (p, 6>(p)) be an element of the set /C. We have: 



9{p) = y/(r + 7 (y)^(p))2 + pV - T , (11) 

with 

H(n\ = . 

i-i{y)%h(r) 
>. n 

8vy 

k 2 v „,v 2 



2 

For y sufficiently large, such that j(y)u = f-^ — > 0, 



8(p) = yfi*Tjw _ r + - T =—H{p)'— + O(-), 

v/r 2 + p 2 v 2 8vy y 2 



and, since H(p) = 4:TcvI (p) + O(^), we obtain the ratio: 



()(/>) \/T 2 + p 2 V — T T , K "1> ^ V 



,2„,2 — a- 1- TC'v V 2 



P P v / r 2 + p 2 v 2 2yp y 2 p 
Therefore, when p — > 0, 

9(p) PV 2 TT 3 U ^. V 2 Up 2 , 

p 2r 2yp |/ 2 p y 



The sum % — h ^ is minimized when p = -, /— , and its minimum is ttv . —. 

It 2yp r v A/ y ' \J V T 

3 

As a result, the ratio ^ is minimized with value ixv^J~-^-~ + ^((y) 2 )' wn i cn corresponds to the 
propagation speed bound. 
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F. Proof of Corollary \2\ 

Again, we take w. 1. o. g., R = 1 and G — 1, and we consider the kernel set (p, 9{p)). From (fTTj) . when 



0(p) = >/(7(!/)^(p)) 2 + A 2 + 0(t). 

We obtain the ratio: 

gM = /(tW^MP + ~ +0(I) 
p V p p 



In this case, y <a<yV ^ v ^ p ^ 2 + t> 2 is minimized when the quantity J(p) = iM^MM [ s a i so minimized. 
We take 7(1/) = J^, since this is an upper bound for any value of the parameters. Thus, using (TT2 
when J 0, 

8t>yp 8yp 



Therefore, the minimum of J(p) is ^ min p ( ; ^^) + O(^), attained for p = 1.608 . . ., and we have the 
propagation speed upper bound: v + 0(^ +r). 



